Cooperative electron-phonon interaction in molecular chains 
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Using a controlled analytic treatment, we derive a model that generically describes cooperative 
strong electron-phonon interaction (EPI) in one-band and two-band Jahn- Teller (JT) systems. The 
model involves a next-nearest-neighbor hopping and a nearest-neighbor repulsion. Rings with odd 
number of sites (o-rings) belong to a different universality class compared to rings with even number 
(e-rings). The e-rings, upon tuning repulsion, undergo a dramatic discontinuous transition to a 
conducting commensurate charge density wave (CDW) state with a period independent of filling. 



PACS numbers; 71.38.-k, 71.45.Lr, 71.38.Ht, 75.47.Lx 

The last few decades have witnessed numerous stud- 
ies to fathom the tapestry of exotic orbital, charge, and 
spin orderings in JT transition metal oxides (TMO) such 
as the manganites, cuprates, nickelates, etc. To 
model the emergent ordering in these complex TMOs, 
one needs, as building blocks, effective Hamiltonians for 
various interactions. Except for the cooperative JT inter- 
action, effective Hamiltonians, that reasonably mimic the 
physics, have been derived for all other interactions. For 
instance, double exchange model approximates infinite 
Hund's coupling, Gutzwiller approximation or dynami- 
cal mean-field theory model Hubbard on-site coulombic 
interaction, superexchange describes localized spin inter- 
action at strong on-site repulsion, etc. Controlled math- 
ematical modeling of cooperative JT quantum systems 
that goes beyond modeling localized carriers 2] has re- 
mained elusive (at least to our knowledge). In fact, a 
controlled analytic treatment of the many-polaron effects 
in single-band Holstein model |3;] (which is a simpler non- 
cooperative EPI system) has been reported only recently 
[j]. However, definite progress has been made in numer- 
ically treating JT systems 

Owing to its cooperative nature, the JT interaction 
leads to non-local distortion effects which can change 
the very nature of long range order. The situation 
is further complicated by the ubiquitous strong EPI 
for which there is compelling evidence in manganites 
(from extended X-ray absorption fine structure [51] and 
pulsed neutron diffraction [6] measurements) , in cuprates 
(through angle-resolved photoemission spectroscopy 0]), 
and in nickelates (based on neutron scattering [§| ) . While 
a weak interaction is amenable to a Migdal-type of per- 
turbative treatment, the strong interaction (even for a 
single-band system) necessitates a non-perturbative ap- 
proach 

In this letter, based on our work 0] on the Holstein 
model [Fig. H^a)], we derive the effective Hamiltonian for 
cooperative strong EPI in one-band system [Fig. [Ijb)] 
and two-band (from d^-i and d^2_y2 orbitals) JT chain 
[Fig. [He)]. Upon inclusion of cooperative effects in the 
strong EPI, the system changes its dominant transport 
mechanism from one of nearest-neighbor hopping to that 
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FIG. 1: Molecular chains with dz2 orbital hopping sites (filled 
circles) and oxygen sites (empty circles) in (a) Holstein model, 
(b) one-band cooperative EPI system, and (c) two-band co- 
operative JT system with Qi, Q2, and Qa octahedral modes. 



of next-nearest-neighbor hopping. Below half-filling and 
at strong EPI in e-rings, while the systems without co- 
operative effects remain disordered Luttinger liquids, our 
cooperative EPI one-band and two-band JT systems pro- 
duce CDW states. Furthermore, like the manganites, our 
JT system too manifests absence of electron-hole symme- 
try only upon including cooperative effects. 

To bring out the essential physics, we begin with a 
one-dimensional (ID) toy-system of electrons hopping in 
a one-band system of d^^ orbitals which are coupled to 
the oxygens in between as shown in Fig. [Ijb) [9|. The 



Hamiltonian is expressed as H = Ht 
the hopping term Ht is given by 
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ff, = -<^(c]c,+i+H.c.), 



Hi where 
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with Cj being the destruction operator of an electron in 
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a d-'a orbital (at site j) and the EPI term Hep by 



with Hi = c\ci, Qi = Ui — Ui-i representing the expansion 
of the oxygens around the d],2 orbital, and the right-hand- 
side (RHS) oxygen displacement = (aj +ai)/^/2MLJo. 
Furthermore, the lattice term Hi representing simple har- 
monic oscillators is of the form 

i j j 

The main difference between the Holstein model and the 
above cooperative Hamiltonian is that in the Holstein 
model electrons at different sites are coupled to different 
on-site molecular distortions whereas in the present toy- 
system the electrons on adjacent sites are coupled to the 
displacement of the same in-between oxygen. Thus in 
our toy-system, to produce an effective polaronic Hamil- 
tonian, we need to devise a modification of the usual 
Lang-Firsov transformation 10] so as to take into ac- 
count the cooperative nature of the distortions. To meet 
this end we used the following canonical transformation 
H = exjp{S)H exp(—S) where S now contains the differ- 
ence in densities on adjacent sites 
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Then, one obtains H = Hq + Hi where 

Ha = 



-te 



H.C.), 



(5) 



and 



Hi=J2 Hi, = -te-'^' 5][c]c,+i{TftT^- - 1} + H.c], 



with 7j! = eKp[±g{2aj 



flj+i)]. On account of 



the cooperative nature of the EPI we obtain an additional 
term 2g'^ujQ njUj^i involving nearest-neighbor repul- 
sion in Hq and the perturbation Hi now involves phonons 
at three sites as opposed to phonons at only two sites 
as in the non-cooperative case. We consider the case 
texp[—3g'^] << ojQ and perform second-order perturba- 
tion theory similar to that in Ref. 0- For large g^, after 
tedious algebra, we obtain the following effective pola- 
ronic Hamiltonian \XU- 

t^ 
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(c]cj+i + H.c.) 



+i(l - nj) 



t^e-^ 



[ct i(l-2n,)c,+i-HH.c.]. (6) 



Notice that the coefficient of the nearest-neighbor hop- 
ping is significantly smaller than the next-nearest- 
neighbor hopping for large g^ and not-too-small i/wg! 
This is a key feature resulting from cooperative effects. 

The above effective Hamiltonian may be contrasted 
with the following Hamiltonian H^ff for the case 
where there is no cooperative EPI [i.e., Hf.p — 
~^/2gu:o n,{al + a,)] [see Ref. i and Fig. [^a)]: 



Heff - 



-te-29'^(c]cj+i +H.C.) 



2g'^uJo 



^[c]_i(l-2n,)c,+i+H.c.].(7) 



Although our results in Eq. ([7]) are valid for te^^^^ << 
ujq and g'^ >> 1, we provide an explanation of the results 
in the more restrictive but physically understandable 
adiabatic (t >> loq) and small polaronic {g^uiQ >> t) 
regime as follows. In Eq. (O, the coefficient of the 
^^ nj-|-i(l — Uj) term can be understood as resulting 
from an adiabatic process when an electron at site j -\- 1 
hops to a neighboring site j and back, but the lattice 
has no time to distort (relax) locally at site j (j -I- 1) 
and thus yields the second-order perturbation energy 
— ^^/(energy change). On the other hand the coeffi- 
cient of ~ 2"'j)cj+i -I- H.c] results when the 
intermediate site j does not distort/relax in an adia- 
batic hopping and thus yields texp[—2g^] x 2pTJ^ where 
iexp[— 2f/^] is due to distortions at both initial and final 
sites j — 1 and j -I- 1. In the above non-cooperative case, 
the nearest-neighbor hopping dominates over the next- 
nearest- neighbor hopping in the small polaron limit. 

Using the above logic we see that the higher order 
terms in perturbation theory, for both cooperative and 
non-cooperative cases, are dominated by the process 
where an electron hops back and forth between the same 
two sites. The dominant term to fcth order is approxi- 
mately given for even k by 



^nj+i(l 



''3h 



while for odd k by 



te- 



t 



fc-i 



5](4c,+i+H.c.), 
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(9) 



where 7 is 2 for the non-cooperative case and 3 for the 
cooperative one. Since each term in the perturbation 
theory should be smaller than wq, we see that the small 
parameter in our perturbation is t/{gujo) [l3]. 

Here, a few observations are in order. Firstly, the coop- 
erative effects, unlike in the Holstein model's case, raise 
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the potential of the site next to an occupied site and 
thus makes it unfavorable for hopping. Thus, in Eq. 
([S]) as compared to Eq. ([7]), the exponent is larger for 
the nearest-neighbor hopping and also the denominators 
of the coefRcients are similarly larger for the hopping- 
generated nearest-neighbor interaction and for the next- 
nearest-neighbor hopping. Next, Lau et al. 0] obtain 
the same energy expression for a single polaron as that 
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given by Eq. ([6]) (when nj ~ 0). Lastly, in Ref. 
the authors explain the ferromagnetic insulating behavior 
in low-doped manganites by using the non-cooperative 
hopping-generated nearest-neighbor interaction [i.e., sec- 
ond term on RHS of Eq. ([7|)] after modifying the hopping 
term for double-exchange effects. From Eq. ([6|), we see 
that cooperative phenomenon must be taken into account 
as it reduces the ferromagnetism generating interaction 
strength by a factor of 1.5. 

We will now analyze the effective polaron Hamiltonian 
given by Eq. We first note that, for large values of 
g^, the essential physics is captured by 

Heff = -7^E(4-i(l - 2n,)c,+i + H.c.) 



(10) 



with T/V << 1. However, owing to its novelty, we will 
study the model given by Eq. pO|l for arbitrary value 
of T/V. Next we note that for T/V « 1 in e-rings, 
the system always has alternate sites (one sub-lattice) 
occupied for less than half-filling and above half-filling 
the other sub-lattice gets filled. This is due to the fact 
that at large repulsion, as we fill up the lattice with elec- 
trons one after the other, each new electron added to the 
system has more number of sites to hop to on the same 
sub-lattice occupied by the previous electrons. As for the 
other extreme situation V = 0, for even number of elec- 
trons on e-rings, the model has both sub-lattices equally 
occupied. 

Using modified Lanczos [3] , we calculated the ground 
state energy using periodic (antiperiodic) boundary con- 
ditions for odd (even) number of fermions. First for e- 
rings with N sites, the ground state energy has a slope 
discontinuity, with the energy increasing up to a criti- 
cal value, after which it is constant (see Fig. We 
will now show clearly that as the interaction strength 
increases, at a critical value of V/T, Ising Zi symme- 
try (i.e., both sub-lattices being equally populated) is bro- 
ken and only a single sub-lattice is occupied for e-rings. 
To this end, we calculate the density-density correlation 
function W{1) = 'i/Nj^jiinj - 0.5){nj+i - 0.5)) and the 
static structure factor S{k) = J^i exp{ikl)W (1) near crit- 
ical value. Now, S'(7r) = [E/cve„ " E/odd]^(0 with 



^(^) = N^^^- - 0-25Af]' + [No ~ 0.257V]2}, (11) 
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FIG. 2: Plots of Energy (E) versus interaction strength (V) 
for rings with A*' sites, Np electrons, and hopping T = 1. 



and 



Zodd 



N 



[Ne - 0.25N][No - 0.25N], (12) 



where iVg(o) is the number of particles on even (odd) sites. 
Then, when Z2 symmetry is respected, for even number 
of particles Np = 2N(. = 2No, we have S{tt) = and for 
odd value of Np we have ^(Tr) — A/N. When only one 
sub-lattice is occupied, on taking N^ = Np and No = 
(without loss of generality), we get ^(Tr) — ANp/N 
and W{lodd) = 1 - 4Np/N. We find that at a critical 
interaction strength, as shown in Fig. [31 the following 
dramatic changes occur: (i) the structure factor S{'k) 
jumps from near to almost its maximum value AN^/N; 
(ii) W{lodd) also jumps close to its large V value of 1 — 
ANp/N; and (iii) W^(/even) (for I ^ 0) too jumps and 
its final value at half- filling is 1. For a fixed T and N, 
the critical value Vf^ of V increases monotonically as Np 
decreases. For N = 16, Np = 2, and T = 1, we get 
Vq^ « 4. From finite size scaling for half-filling, using 
V^ -V^ (X 1/N'^ and system size N < 20, we obtain 
V(f « 2.83. 

We see from the above analysis that, at a critical repul- 
sion, the system undergoes a discontinuous transition to 
a conducting commensurate CI? state away from half- 
filling while at half-filling one obtains a Mott transition. 
Usually commensurate CDWs are insulating (see Ref. 
[lih whereas our model surprisingly predicts a conduct- 
ing commensurate CDW. Furthermore, quite unlike the 
Peierls transition, the period of the CDW is independent 
of density] Such density independent charge ordering has 
indeed been observed in manganites (see Fig. 2 in Ref. 
[l6[ ) . Our toy-system is different from the non-cooperative 
EPI Holstein model in that the latter, at strong cou- 
pling and finite t/coo, can be mapped onto an anisotropic 
Heisenberg model and is thus always a Luttinger liquid 
away from half- filling (see Ref. 0). 

We will now discuss o-rings for number of electrons 
less than half the number of sites. The case when the 
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FIG. 3: Plot of the density-density correlation function W{1) 
and structure factor S{k) for a 16-site ring with Np electrons, 
interaction strength V, and hopping T = 1. 



number of electrons is greater can be deduced invoking 
electron-hole symmetry. Since the system has all sites 
connected through next-nearest-neighbor hopping (much 
like a Moebius strip), there do not exist 2 sub-lattices. 
Consequently, as the interaction increases, the system's 
energy increases monotonically and smoothly (without 
any slope discontinuities) and there is no phase transition 
[see Fig. [5]. We find that the system's energy increases 
linearly at small strengths as can be expected from a 
mean- filed analysis. At large values, the system's energy 
approaches asymptotically the energy of the e-ring with 
one less site but with the same number of electrons. This 
is because at large V, hard core repulsion on adjacent 
sites blocks electrons similar to Pauli's exclusion principle 
on the same site. 

Lastly, for the ID JT system shown in Fig. [ijc), we 
first note that there is a strong on-site inter-orbital repul- 
sion given by Uj2j '^j^'j where U ^ oo and n'j = djdj 
with d being the destruction operator of a localized elec- 
tron in a dx2_y2 orbital. Next, like the one-band case [see 
Eq. ([6])], the cooperative JT effect too produces strong 
nearest-neighbor replusion oc g^uJa^^tijUj+i. Again 
as before, there is a virtual-hopping generated nearest- 
neighbor repulsion term —G^- g^j^-y nj+s(X — nj)(l — n'j) 
with the additional (1 — n'^) factor prohibiting hopping 



of a d^-2 electron to a neighborging j site if it is occupied 
by a dx2_y2 electron. Here too G oc /{g^ujo). Then, the 
Hamiltonian at large EPI (whose exact form is obtained 
after quite tedious algebra [Hj) is essentially given by 

^eff = + Un^n'^] -G - "^■) 

-r^[(l - n'^)c]_,{l - 2n,)c,+i + H.c], (13) 



with G » T. As in Eq. PH)) . here too V is an arbitrary 
variable. In Eq. (jl3p . in the next-nearest-neighbor hop- 
ping term, the (1 — n'j) factor does not allow hopping via 
the i site if it is occupied by a dx2-y2 electron. 

For fiUings up to half (i.e., 0.5 electron/site), the 
electrons populate only the dz2 orbitals (i.e., n'^ = 0) 
and the behavior of the above model of Eq. (fT3|) is 
identical to that of the effective model given by Eq. 
PU)) and is depicted by Fig. [5] when E is replaced by 
E -\- 2NpG. Our results agree with the experimental ob- 
servations that C-type ferromagnetic-chain manganites, 
like Ndi^^Sr^MnOs at 0.63 < x < 0.8 are d^2 

polarized. For fillings between half and unity and when 
V >> G {as in a real cooperative EPI system), for e-rings 
[o-rings] N/2 [{N — l)/2] electrons populate orbitals 
with the rest occupying dx2_y2 orbitals; thus system lacks 
electron-hole symmetry. 

On ignoring cooperative effects in our JT sytem of Fig. 
[Ijc), in Eq. ([T3|) we set V = 2G and add a nearest- 
neighbor hopping -To J2j,s=±i (^l^j+s with G >> Tq >> 
T [as in Eq. ([7])]. Then, importantly, only dz2 orbitals 
are occupied at any filling (up to unity) and the sys- 
tem exhibits electron-hole symmetry. Furthermore, since 
n'j = 0, the system is always a Luttinger liquid at non- 
half-fiUing 4]. 

In summary, our study shows that cooperative EPI 
produces strikingly different physics compared to the 
non-cooperative situation. Our derived model captures 
the essential feature of cooperative strong EPI in all di- 
mensions (and in even non-cubic geometries), i.e., the 
dominant transport is due to next-nearest-neighbor hop- 
ping. Our controlled analysis of cooperative JT interac- 
tion is currently being extended to higher dimensions. 
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